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Based upon the rate equations for the photon distribution function obtained in the previous paper, 
we study the inverse Compton scattering process for high-energy nonthermal electrons. Assuming 
the power-law electron distribution, we find a scaling law in the probability distribution function 
Pi(s), where the peak height and peak position depend only on the power index parameter. We 
solved the rate equation analytically. It is found that the spectral intensity function also has the 
scaling law, where the peak height and peak position depend only on the power index parameter. The 
present study will be particularly important to the analysis of the X-ray and gamma-ray emission 
models from various astrophysical objects such as radio galaxies and supernova remnants. 
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I. INTRODUCTION 

The inverse Compton scattering is one of the most fun- 
damental reactions which have variety of applications to 
astrophysics and cosmology. They are, for example, the 
Sunyaev-Zeldovich (SZ) effects[l| for clusters of galax- 
ies (CG), cosmic-ray emission from radio galaxies [2| and 
clusters of galaxies |3j, and radio to gamma-ray emission 
from supernova remnants [3, Therefore, theoretical 
studies on the inverse Compton scattering have been 
done quite extensively for the last forty years, start- 
ing from the works by Jones @, and Blumenthal and 
Gould to the recent works, for example, by Fargion[l], 
Colafrancesco(9l IToj . and Petruk(TT|. 

In particular, remarkable progress has been made in 
theoretical studies for the SZ effects for CG. Wright [12f 
and Rcphacli 13] calculated the photon frequency redis- 
tribution function in the electron rest frame, which is 
called as the radiative transfer method. On the other 
hand, Challinor and Lasenbvflij and Itoh, Kohyama, and 
NozawaflE) solved the relativistically covariant Boltz- 
mann collisional equation for the photon distribution 
function, which is called the covariant formalism. Al- 
though the two are very different approaches, the ob- 
tained results for the SZ effect agreed extremely well. 
This has been a longstanding puzzle in the field of the 
relativistic study of the SZ effect for the last ten years. 
Very recently, however, Nozawa and Kohyama[16j (de- 
noted NK hereafter) showed that the two formalisms were 
indeed mathematically equivalent in the approximation 
of the Thomson limit. This explained the reason why the 
two different approaches produced same results for the SZ 
effect even in the relativistic energies for electrons. 



In the present paper, we extend the formalism obtained 
by NK to the case of high-energy electrons. This exten- 
sion will be particularly interesting for the analysis of 
X-ray and gamma-ray emissions, for example, from ra- 
dio galaxies[2j and supernova remnants 0, where the 
inverse Compton scattering of the CMB photons off non- 
thermal high-energy electrons plays an essential role. In 
the present approach, we push analytic techniques as 
much as possible in order to obtain analytic solutions. In 
contrast to the direct numerical calculation, the present 
approach will have an advantage that one may reveal es- 
sential physics properties behind the numerical results. 
In the present paper, under a specific condition for the 
electron distribution which is typically realized, we will 
show that a universal scaling law is valid for the spectral 
intensity function. 
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The present paper is organized as follows: Starting 
from the rate equations derived in the NK paper, we 
derive in Sec. II the analytic expressions for the redis- 
tribution function P(s, 7) and probability distribution 
function P\(s). Assuming the power-law electron dis- 
tribution, we show that -Pi(s) has a scaling law, where 
the peak height and peak position depend only on the 
power index parameter. We calculated the rate equation 
and obtained the analytic expression for the spectral in- 
tensity function dI(X)/dr. We show that dI(X)/dr also 
has the scaling law, where the peak height and peak po- 
sition depend only on the power index parameter. In 
Sec. Ill we apply the scaling law to the observation of 
the spectral intensity in the X-ray and gamma-ray en- 
ergy regions. Finally, concluding remarks are given in 
Sec. IV. 
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II. HIGH-ENERGY INVERSE COMPTON 
SCATTERING 

A. Rate Equations in Thomson Approximation 

In the NK paper, it was shown that the covariant 
formalism [l5| and radiative transfer method|l2| were 
mathematically equivalent in the following (Thomson 
approximation: 



follows: 



7— < 1 

TO 



1 



Vw 2 "' 



(i) 

(2) 



where w is the photon energy, 7 is the Lorentz factor, and 
j3 and to are the velocity and rest mass of the electron, 
respectively. Throughout this paper, we use the natu- 
ral unit h = c = 1, unless otherwise stated explicitly. 
For the cosmic microwave background (CMB) photons 
{k B T C MB = 2.348 x lCT 4 eV), u) < 5 x l(T 3 eV is well 
satisfied. Then w/m<lx 10 -8 , which implies 7 <C 10 s . 
Therefore as far as the CMB photons are concerned, 
Eq. ([1]) is fully valid from non-relativistic electrons to 
extreme-relativistic electrons of the order of TeV region. 

The rate equations for the photon distribution function 
n(x) and spectral intensity function I(x) were derived in 
the NK paper under the assumption of Eq. ([1]). Here, 
x = uj /ksTcMB is the photon energy in units of the 
thermal energy of the CMB. We recall the results here 
to make the present paper more self-contained. They are 
given as 

follows [Mill: 



dn(x) 

dr 
dl(x) 

dr 



dsPi(s) [n(e s x) - n(x)] , 



(3) 



dsP 1 (s)[l(e- s x)-I(x)] , (4) 
t = n e a T t , (5) 



where I(x) — I x 3 n(x), I — (k B TcMB) 3 /27r 2 , n e is the 
electron number density, er<r is the Thomson scattering 
cross section. In Eqs. ([3]) and d4|, Pi(s) is the probability 
distribution function for the photon of a frequency shift 
s, which is defined by e s = x' /x, 



Px(s) = 



Pmin 

e s 



dpp 2 f Pe {E)P{s,P), 

'M2(s) 



2/?7 4 



/(Mo,Mo) = o 



1 



(6) 

/(Mo, Mo) . ( 7 ) 
(8) 



where p e (E) is the electron distribution function of a mo- 
mentum p which is normalized by J °° dpp 2 p e (E)/m 3 = 1. 
Variables appearing in Eqs. © - © are summarized as 



/3 mm = (l-e-l s ')/(i + e" |s| ) 
M = [l-e s (l-/3/x )]/£, 



(9) 
(10) 



/ \ J -1 for s < , . 

mS> ~ \ [1 -e- s {l +ff)]/0 fors>0 ' { ' 



by 



The total probabilities for P(s,j3) and -Pi(s) are given 



dsP(s,j3) = 1, 

A3 

/OO 
dsP^s) = 1, 
-OO 



(13) 
(14) 



where 



Xp = In 



1 + /? 
1-/3 



(15) 



It should be noted that the following useful relations 



P(s,f3)e- 3s =P(-s,0) : 
Pi(s)e- 3s = 



(16) 
(17) 



are valid. 



B. P(s,/3) for Extreme-Relativistic Electrons 

In this section, we derive the analytic expression of 
the frequency redistribution function P(s, j3) for extreme- 
relativistic electrons. In Eq. 0, the integral of [1q can 
be done analytically. One obtains as follows: for s < 0, 



P(s,/3) 



[-Ci^-CaCSJe' + Cs^) 



,2.s 



32/3 2 7 4 

+C 4 (/3)(A^+ s )( e s + e 2s ) + C 1 (/3) e 3s ] , (18) 



and for s > 0, 



[CiG8) + C 3 G8)e'-C a (/3)< 



2.1 



32/3 2 7 4 

+C 4 (/3)(A, 3 - S )(e s +e 2s )-C 1 (/3) e 3s ] , (19) 



where the coefficients are 
1 



Cxifi) 



/? 4 7 2 ' 
1 



(20) 



W) = pA{1 - (4/3 4 ~(3 3 ~ 13/3 2 - 3(3 + 9) , (21) 
C 3 (/3) = /?) + P 3 ~ + 3/3 + 9) , (22) 

C 4 (/3) = ^(/? 2 -3) . (23) 
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Note that Eqs. (118)) and (1191) satisfy the relation of 
Eq. (fTOj) . It should be also noted that Eqs. (fTS|) and 
(fT9f agree with Eqs. (23a) and (23b) of Fargion et al.@, 
respectively. 

Now let us consider the case for electrons of extreme- 
rclativistic energies E (= 7771c 2 ) 3> mc 2 . Thus, 7 > 1 
and P ~ 1 are assumed. In this approximation, Eqs. ([l£ 
and (|19p are written as follows: for s < 0, 



32 7 4 



+ 2e s + 8 7 2 e 2s 



and for s > 0, 



P(*,7) 



32 7 4 



-4(A 7 + s) (e s + e 2s ) + 



+ 87V + 2e 2s 



„3s 



(24) 



-4(A 7 -s) (e s + e 2 *)-- 



,3s 



A 7 = 21n(2 7 ) , 



(25) 
(26) 



where the expression P(s, 7 ) was used instead of P(s, 0). 
Equations ([24]) and |25|) can be integrated analytically. 
One obtains as follows: for s < 0, 



/-> p( " ) = (?) 



and for s > 0, 



rfsP(s, 7 ) = 1 + ( — 



(27) 



(28) 



In Eq. (|28|) . the terms contributing to O (l/7 2 ) are the 
higher-order terms. Therefore one can eliminate the cor- 
responding terms from Eq. (|25p , which gives the new ex- 
pression for s > 0. Then the new expression for s < 
is obtained with the relation of Eq. (fT6|) . Therefore, the 
total probability 



J + \sP{s, 1 ) = l + o( y ^ 



(29) 



is satisfied for P(s, 7 ). The explicit forms are as follows: 
for s < 0, 



Let us now compare the present results with the litera- 
ture. It is straightforward to show that Eqs. ([501) and (|3"Tj) 
are equivalent to Eqs. (38) and (40) of Jones[6j, respec- 
tively. We show the equivalence between the present for- 
malism and Jones's formalism in Appendix A. It should 
be also mentioned that Eqs. (24a) and (24b) of Fargion 
et al.Q differ from our Eqs. fl3U} and ([5T|). The difference 
comes from 0(l/ 7 2 ) terms as mentioned in their paper. 

Before closing this subsection, it should be also noted 
the following: In the present formalism, the cases s > 
and s < correspond to the Compton scattering and in- 
verse Compton scattering, respectively. This is because 
of the definition x — e~ s x' , where x' and x are the en- 
ergies (in units of }zbTcmb) of initial and final photons, 
respectively. Equations (|2"7|) and ([2"5]) suggest that prob- 
ability distribution for the CMB photon scattering by 
high-energy electrons is dominated by the Compton scat- 
tering process instead of the inverse Compton scattering 
process. 



C. Scaling Law of Pi(s) for Nonthermal Electrons 

In order to proceed calculation for practical applica- 
tions, let us specify the electron distribution function. 
High-energy electrons in the supernova remnants and ac- 
tive galactic nuclei, for example, are most likely nonther- 
mal. It is standard to describe the nonthermal distribu- 
tion in terms of the power-law distribution function of 
three parameters: 



elsewhere 



(32) 



where 7 is the Lorentz factor and A 7 is the normalization 
constant. In Eq. Q32p. a is the power index parameter, 
7mm and 7maa; are parameters of minimum and maxi- 
mum values for 7, respectively. Then, Eq. © can be 
reexpressed as follows: for s < 0, 



Px(s)= / d 7 P e (7)^, 7 ), (33) 

Jmax(7 mi „,e- s /!/2) 

where P(s, 7) is given by Eq. (|5U|) . and for s > 0, 



and for s > 0, 



32 7 4 



1 



+ 2e s + 87 2 e 2s 



32 7 4 



-4 (A 7 + s) e 

87V + 2e 2s 

-4 (A 7 - s) e 2s - • 



„3s 



(30) 



(31) 



Pi (»)= / d 1Ve {i)P{s,i) , (34) 

Jmax(7 mi „,e s / 2 /2) 

where P(s, 7) is given by Eq. pip . In deriving Eqs. (j33p 
and (pH)) , /3 « 1 was assumed, and the phase space factor 
7 2 was absorbed, for simplicity, by the power index a in 
Veil)- 

In the case of the power- law distribution of Eq. (]32l) . 
equations (|3"3"|) and (j34D can be integrated analytically. 
The explicit forms are given as follows: for —2 ln2"/ rnax < 



4 



s < -21n27 mi „, 

62 I cr + 5 



Jmax 



(7 + 3 



£7- 1 
cr + 3 

1 /cr - 1 



q(7+3 (ct+3)s/2 
I 

- 2s-41n27„ 



cr + 1 

for -21n27 mi „ < s < 0, 
1 



7mtx V CT + 3 

2<T+l e ( ( T+l) S /2 



1 



Jmax 



(35) 



_3_ 

32" 
2 

(7 + 3 



^+3 
/ mm 



/T + 3 
max 

l 



f 1 


5 


^ imin 


(7 


- 1 


£7 


+ 3 


(7 


- 1 


(7 


+ 3 



1 



-2s - 41n27„ 



- 2s-41n27„ 



. 2.s 



a + 1 \ 7 TO "| n 7moa 
for < s < 21n27 mi „, 

1 / I 



(36) 



+ 



3 

32" 
2 

cr + 3 



1 



• ^min ^max 
1 /cr - 1 



„3s 



2s - 4 In 27„ 



7 CT + 3 I cr + 3 

/ mm V 
1 /cr - 1 



,<r+3 I _ , o 
max \ w 1 u 

1 I 



2s - 41n27„ 



cr + 1 l7 CT+1 7mti 

\ hmn jmax 

and for 21n27 min < s < 2h\2^ max , 
3 



(37) 



cr + 3 



CT + 5 
cr - 1 



1 2 <T+5 e -(°'+ 5 ) ;5 / 2 ' 



Jrnax 



Jmax 



2<H-3g-(er+3)s/2 

2s - 41n27„ 



cr + 3 

1 /cr - 1 
Cr + 3 



_ 2 s 



+ _8_ / 2CT+ l e -( ff +l)s/2 _ 1 
^+1 V 7mti 



(38) 



It should be noted that the normalization constant is 
given by 



iV 7 = (cr - 1) 7 ; 



<x-l 
min 



(39) 



for the case 7 max — > °o. 

Let us now introduce new functions Pc(s,R) and 
Pj C (s,R) in order to express Eqs. (|33|) - (l58l in unified 
forms, where R = ^min/lmax- Here, C and IC denote 



the Compton scattering and Inverse Compton scatter- 
ing, respectively. First, we define Pc(s,R) as follows: 
for -21n27 mi „ < s < 0, 



Po(s,R)=3 



cr - 1 

1 - R°~ 
1 

+ 

cr + 3 

_ R *+3 

1 



(7 + 5 



cr - 1 
cr + 3 
cr- 1 
cr + 3 " 



-2s 

2s + 4 In R ) 



"1 - e s 



Cr + 1 

and for < s < 21n(7 moa /7 mi „), 



(40) 



P c (s,R) = 3 



cr- 1 



1-R° 



a + i> 



e -( CT +5)s/2 _ ^+5} e 3s 



1 



cr + 3 



cr - 1 
cr + 3' 



-(<t+3)s/2 



-i^+s f £ — - + 2s + 41ni2 

cr + 3 > 



„2s 



+ _l_( e -(-+l) S /2_ i?CT +l) e s| _ 

Similarly, P' w (s,R) is for -2 ln(7 maa; /7 min ) < s < 



(41) 



1 



cr + 3 
1 



,(<t+3)s/2 



cr + 1 

and for < s < 2 In 2^ m i n , 
cr - 1 



cr - 1 
cr + 3 L 

Cr- 1 \ 

— — -2s + 41ni? 

cr + 3 / 

e (<r+l)s/2 _ R a+A £ 2 



(42) 



p; c (s,i?) = 3 



1 - R 
1 

"cr + 3 
_ R *+3 



er-1 



a + t> 



- (i - ir+ 5 ) 



(7 - 1 
cr + 3 
cr - 1 
cr + 3 



- 2s 

2s + 4 In .R | 



It is straightforward to show that 

Pc(s,i?)e- 3s -P; c (-s, J R) 



(43) 
(44) 



is satisfied by Eqs. (|5D ]) -(|i5 |l . 

Comparing Eqs. (p ]) -(|38 ]l with Eqs. (|gQ )) -(|i3 ]l , the 
probability distribution function Pi(s) is described as fol- 
lows: 



Pi(s) 



Pic (s + 2 In 2 7mm ,R) for s < 
P c (s- 2 In 2-y min ,R) fors>0 



(45) 



5 



where 



Pic(s,R) 



64 7 ; 



6 

rain 



Pic(s,R). 



(46) 



Let us now consider the case R = Jmin/lmax "C 1. We 
fix j max = 10 throughout the paper. In Fig. 1(a), we 
plot -Pi(s) as a function of s for a typical value a = 2.5. 
The solid curve, dash-dotted curve, dashed curve, and 
dotted curve correspond to 7 m i„ = 10, 10 2 , 10 3 , and 
10 4 , respectively. It can be seen that the height of Pi(s) 
is independent of 7m i„. In Fig. 1(b), we plot the same 
curves as a function of new variable sc which is defined 
by 



sc 



21n27„ 



(47) 



In Fig. 1(b) the four curves are totally indistinguishable, 
which exhibits a scaling law for Pi(s). The reason for this 
scaling law is as below. For large 7 m i„ 3> 1, as shown by 
Figs. 1(a), 1(b), and Eqs. g3]) and (pfg]h the probability 
distribution function Pi(s) is dominated by Pc(sc,0), 
i.e. by the Compton scattering process. 

Before closing this subsection, we study the a- 
dependences on the peak position s pea k and peak height 
Pi(speak)- As shown in Figs. 1(a) and 1(b), the j m in- 
dependence of Pi(s) is described by Eq. (|4"?| . namely, s 
= sc , + 21n27 m i„. Therefore, we define the peak position 
by 



Speak = s(a) + 21n27„ 



(48) 



where s(a) depends only on a. The peak position is 
calculated by solving the equation 



dPi(s) 



o. 



(49) 



The analytic expressions for s(a) in the first-order and 
third-order approximations are given as follows: 



sist(o-) : 
S3rd(cr) 



((T-l)((7 2 +4(7+ll) 



5a 3 + 23a 2 
1 



(50) 



2(4o- 2 + 21cr + 29) 



51cr+ 17 ' 

3cr 2 + 14c7 + 19 (51) 



y/(a + I) 2 A 3 + B 2 + B 
(J + 1 / 
' y/(cr + l) 2 A 3 + B 2 -B* 

(7+1 



1/3 



1/3' 



. (52) 



A = 7<r i - 
B = 3cr 7 



- 64a 3 + 254a 2 + 520a + 451 , 

- 21cr 6 - 582cr 5 - 4378cr 4 - 18589ct 3 
^48333cr 2 - 70688cr - 44036 . 



(53) 



(54) 



We also solved Eq. (|4"9")l numerically and obtained the 
numerical solution s num (a). In Figs. 2(a) and 2(b), we 




FIG. 1: Plotting of Pi (s) and Pi (sc) for a = 2.5. Figures 1(a) 
and 1(b) are Pi(s) and Pi (so), respectively. The solid curve, 
dash-dotted curve, dashed curve, and dotted curve correspond 
to jmin = 10, 10 2 , 10 3 , and 10 4 , respectively. 



plot s(a) and Pi(s pea k), respectively. The dashed curve, 
dash-dotted curve, and solid curve correspond to si s t{a), 
S3rd(c) and s num (a), respectively. In Fig. 2(b), the solid 
curve and dash-dotted curve are indistinguishable. It can 
be seen from Figs. 2(a) and 2(b) that the third-order ap- 
proximation is sufficiently accurate for the present pur- 
poses. 



D. Scaling Law for Spectral Intensity Function 

Let us now solve the rate equations of Eqs. ([3]) and 
d3| with the result of Eq. (gSJ) for P%{s). We consider 
the CMB photons for the initial distribution. For the 
inverse Compton scattering by high-energy electrons, we 
are interested in high-energy spectrum such as X-rays 
(~ keV) and gamma-rays (~ MeV). Therefore, one can 
safely assume 



> 1 



(55) 
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FIG. 2: Plotting of s(a) and Pi(s peak ). Figures 2(a) and 2(b) 
are s(a) and Pi(s pea k), respectively. The dashed curve, dash- 
dotted curve, and solid curve correspond to the first-order 
approximation, third-order approximation, and numerical so- 
lution, respectively. 



for scattered photons. For the 7-parameters, we assume 
the same condition used in the scaling law for Pi(s), 
namely, 



1 < 7mm < 7 II 



(56) 



Under these assumptions, Eqs. ((3|) and (|4]) are much sim- 
plified, and can be solved analytically. The derivation is 
straightforward, however, it is lengthy. Therefore, we 
give the derivation in Appendix B in detail, and quote 
the final results here. 

According to Eqs. (IB38|) and (|B40j) , one has the fol- 



lowing results: 
dI{X) 



dr 



3(cr 
1 



1)X S 



x 



+ 



f 3 

6(cr - l)(a 2 +4a 



a a 

T 2 



a- 1 



11) 



dt 1 

Te* - 1 
t \ * 

21n x )x 



1 



(<t + 1)(ct + 3) 2 (ct + 5) ^' 2 
dn(X) 1 1 dI(X) 



x 



dt 



1 t 2 
+ a + lX 2 

t (<r+3)/2 



1 



dT 

X 



64 7 ; 

x 

rnin 



6 

rnin 



I X 3 dT 



,(57) 

(58) 
(59) 



where Iq = (fc^TcMs) 3 /2-7T 2 . It should be empha- 
sized that the function dI{X)/dr depends on 7 m i n only 
through X . Therefore, dI(X)/dr has the scaling law. On 
the other hand, the function dn{X) / dr does not have the 
scaling law because of the factor 1 / 647^ in in Eq. . 
For X>1, Eq. ([57)) is further simplified as follows: 



dl(X) 
dr 



= ^0 

xr 
1 



3(cr-l)(cr 2 +4cr + ll) 

(<j + l){a + 3) 2 (a + 5) 
<y + 5\ Vcr + 5 s 



C 



X 



(<T-l)/2 



m Jo 



dt- 



1 



(60) 
(61) 



where ((z) is the Riemann's zeta function. We call 
Eq. (|60j) the power-law approximation. 

In Fig. 3, we plot dI(X)/dr of Eq. (|57| as a function 
of X for typical cr-values for illustrative purposes. The 
dashed curve, dash-dotted curve, and solid curve corre- 
spond to cr=2.5, 3.5, and 4.5, respectively. The peak po- 
sition and peak height depend only on the power-index 
parameter. It should be noted that dI(X)/dr cx X for 
X < 1, and dI(X)/dr cx X^ a '^l 2 for X > 1. There- 
fore, the slope of the downward curves in Fig. 3 will de- 
termine the cr-value. 

In Figs. 4(a) and 4(b), we plot the peak position and 
peak height of the spectral intensity function as a func- 
tion of <7, respectively. The solid curves correspond to the 
numerical values. The dash-dotted curves are the results 
of analytical fitting functions. They are given by 



Xpeak — 
dI(Xp ea k) 

drl 



1 



a + aiS 1/4 + a 2 S 1/2 

_ 3 E(4 + 6S + E 2 ) 
~ I^o + oxE + ^ + E 3 ) ' 
a- 1. 



(62) 

(63) 
(64) 



The fitting parameters are ao=— 2.18351, a\=h. 37131 
and ct2=— 2.02638 for the peak position, and 6o=2. 60331, 
&i=6.6352 and &2=5.6526 for the peak height. The er- 
rors of the fitting functions in the region 2 < a < 10 are 
less than 0.15% and 0.10% for X pea k and dI(X pea k)/dT, 
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FIG. 3: Plotting of dI(X)/dr. The dashed curve, dash-dotted 
curve, and solid curve correspond to a = 2.5, 3.5, and 4.5, 
respectively. 



respectively. In Fig. 4, two curves are totally indistin- 
guishable. 

Before closing this section, let us compare the present 
result with the literature. In the textbook by Rybicki and 
Lightman[l8j. the expression for the power-law approxi- 
mation is given by Eq. (7.31). The scaling law is hidden 
in the expression of Eq. (7.31). Inserting the explicit 
form of the normalization constant C = n e (a — 1)7^7™ 
in Eq. (7.31), one finally obtains as follows: 



dE 



dVdrdei 



oc I G{a)X- [a - 1)/2 , 



G(a) = 



6(cr- 1)((T 2 + 4(7 + 11) 



Xf 



(<7 + l)(cr + 3) 2 ( ( J + 5) 
<T + 5\ „ /(J + 5" 



(65) 



(66) 



which agrees with Eq. (|60)l . The applicable range of the 
power-law approximation of Eq. (7.31) is given [18j by 



4 7 



2 

min 



El 2 

^ ^ ^7mai ' 



(67) 



where e is a typical energy of initial photon distribution. 
In the case of the CMB photon distribution, one can use 
e = ksTcMB- Therefore, one obtains the condition for 



2.5 
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2 3 
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0.7 



0.6 - 



0.5 



4(b) 



Peak Height dI(X peak )/dr/I 




FIG. 4: Plotting of the peak position and peak height of the 
spectral intensity function as a function of a. The solid curves 
correspond to the numerical values. The dash-dotted curves 
are the results of the analytical fitting functions. 



X as follows: 



R 2 



(68) 



which again agrees with the condition of the present pa- 
per. It is needless to mention that the full expression of 
Eq. (I5T)) has to be used for X < 0(1) as shown in Fig. 3. 



III. ASTROPHYSICAL APPLICATIONS OF 
SCALING LAWS 

A. X-ray region 

In the present section, we show an application of the 
scaling law. Recently, observational studies on the in- 
verse Compton scattering have been done quite exten- 
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5(a) 



5(b) 




5(c) 



5(d) 




Eq. (|63[) . This will serve as an independent check for 
the CT-value. On the other hand, the 7 mm -value is deter- 
mined by measuring the peak position uj pea k in Fig. 5. 
Using the relation of u> pea k with X pea k , namely, 



X 



1 



peak 



4 7mi„ k B T C MB ' 



the 7„ 



-value is determined by 

1 (dpeak 



7mm 



4X„ ea fe ksTcMB 



1/2 



(70) 



(71) 



where X peak is calculated by the RHS of Eq. (Ifj2"j) with 
the measured a- value. It can be seen from Fig. 5 that 
the X-ray observations have sensitivities to 7mm =500 ~ 
3x 10 3 region. 

Before closing this subsection, let us study the appli- 
cability of the power-law approximation used in the lit- 
erature. The condition for the power-law approximation 
X » 1 reads 



(72) 



In the case of j m i n — 1 x 10 3 , for example, one has u> 3> 1 
keV. As shown in Fig. 5(b), the error of the power-law 
approximation is quite large in uj ~ 0(1) keV region. 



0.01 1 — — — 1 0.01 1 — — — 1 

0.1 1 10 100 0.1 1 10 100 

oj (keV) a (keV) 

FIG. 5: Plotting of dI(X)/dr as a function of the photon 
energy cj in X-ray energy regions for a typical value u=2.5. 
Figures 5(a), 5(b), 5(c) and 5(d) correspond to 7 min =500, 
lxlO 3 , 2xl0 3 and 3xl0 3 , respectively. The solid curve cor- 
responds to the full calculation of Eq. (|57[) . The dash-dotted 
curve is the power-law approximation of Eq. (|60[1 . 

sively, for example, X-ray observations from radio galax- 
ies with Chandra[l9ll20j. 

In Fig. 5, we plot dI{X)/dr as a function of the pho- 
ton energy u in X-ray energy region for a typical value 
cr=2.5. Figures 5(a), 5(b), 5(c) and 5(d) correspond to 
7 min =500, lxlO 3 , 2xl0 3 and 3xl0 3 , respectively. The 
solid curve corresponds to the full calculation of Eq. (|57| . 
The dash-dotted curve is the power-law approximation of 
Eq. (|60l) . In Fig. 5, the peak height is independent of the 
7mm -values as pointed in the last section. On the other 
hand, the peak position is shifting toward to high-energy 
side as the j m i n -value increases. 

By measuring the slope of the downward curve in 
Fig. 5, one can determine the er-value, because 

oc ^-iy* (69) 

is valid. One can also determine the cr-value by mea- 
suring the peak height in Fig. 5 with the expression of 



B. gamma-ray region 

With the scaling law for the spectral intensity func- 
tion, one can extend the present formalism to the gamma- 
ray region. In Fig. 6, we plot the same figure as Fig. 5 
for the gamma-ray region. Because of the scaling law, 
the factor vl000(w 31.6) should be simply multiplied to 
the 7 mm -values of the keV region in order to obtain the 
spectral intensity function in the MeV region. Therefore, 
observations in this energy region will have sensitivities 
to 7„ lm =16xl0 3 ~ 95xl0 3 region. Similarly, the factor 
1000 should be multiplied to the 7„ lm -values of the keV 
region in order to obtain the parameter values in the GeV 
region. 

IV. CONCLUDING REMARKS 

In the NK paper [l6|, we derived the frequency redistri- 
bution function P(s, f3) for a frequency shift s and elec- 
tron velocity /3. The form was derived in the Thomson 
approximation, however, it was mathematically equiva- 
lent to that in the covariant formalism [15j. Therefore the 
frequency redistribution function can be applicable from 
nonrelativistic electrons to extreme-relativistic electrons. 

In the present paper, we have extended the formal- 
ism to extreme-relativistic electrons. First, we derived 
the analytic expression for P(s, 7) in the approximation 
7 > 1, It has been found that the present formalism is 
equivalent to Jones's formalism Q. 
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6(a) 



6(b) 




1 10 

cj (MeV) 



FIG. 6: Plotting of dI(X)/dr as a function of the photon 
energy uj in gamma-ray energy regions for a typical value 
a=2.5. Figures 6(a), 6(b), 6(c) and 6(d) correspond to 
7mm=15.8xl0 3 , 31.6xl0 3 , 63.2xl0 3 and 94.9xl0 3 , respec- 
tively. The solid curve corresponds to the full calculation of 
Eq. (1571) . The dash-dotted curve is the power-law approxima- 
tion of Eq. Ml. 



By averaging P(s, 7) over the nonthermal electron dis- 
tribution function, we have calculated the probability dis- 
tribution function -Pi(s). As for the nonthermal distri- 
bution function, we have adopted a standard power-law 
distribution function of three parameters: the power in- 
dex cr, minimum value jmin, and maximum value j ma x of 
the distribution range. For the case j m in 3> lj we have 
found a scaling law in Pi(s), where the peak position 
depends on s — 2ln2j m i n , and the peak height depends 
only on the power index parameter a. 

We have calculated the spectral intensity function. For 
the case of high-energy photons of x ^> 1, we have found 
a scaling law in dl(x)/dr, where the function depends on 
a new variable X = x/(4r/^ nin ). The peak position and 
peak height depend only on the power index parameter 
a. The 7 m i„-dependence of dI(X)/dT is included in the 
variable X. 

We have applied the present formalism to the observa- 



tion of the spectral intensity function in the X-ray and 
gamma-ray energy regions. It has been found that the 
sensitivities of the observation in the X-ray and gamma- 
ray regions are 7mm=500 ~ 3xl0 3 and 7 m i„=16xl0 3 ~ 
95xl0 3 , respectively. 

Finally, we have studied the applicability of the power- 
law approximation used in the literature. In the case of 
Imin = 1 X 10 3 , for example, the error of the power-law 
approximation is quite large in uj ~ 0(1) keV region. 
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Appendix A: Comparison with Jones's Formalism 

The double differential cross sections for cxtrcmc- 
rclativistic electrons are given by Eqs. (38) and (40) of 
Jones's paper (> as follows: 



d 2 N 


2irr 2 c 


dtda 


a-il 2 


d 2 N 




dtda 


2a 17 4 



2q\nq+{l + 2q)(l-q) 



2 In q' 



(Al) 
(A2) 



The variables in Eqs. (IA1[) and (|A2[) are related to the 
variables of the present paper as follows: 
a/ax = e~ s , q = e~ s /47 2 , q' = 47 2 e~ s , a = 9 C mbx, 
where Ocmb = ^bTcmbI^c 1 ■ With these variables, 
Eqs. (jAll) and (|A2| are rewritten as follows: 



d 2 N _ 3a T c 1 
dtdx ~ 327 4 a; 6 



-3s 



7 

■4 (A 7 + s) e" 



8 7 2 e 2s 



d 2 N 



dtdx 327 4 x 



2.s 



87V + 2e 
-4 (A 7 - s) e 25 ' - ^e 3s 



(A3) 



(A4) 



Let us denote the photon distribution function in 
Jones's formalism as nj(a\). Then one has 



nj(ai)dai 



-da\ 



7r 2 (ftc) 3 e ai l 9cMB - 1 

^3 

-x A e As n(e s x)ds , 



(k B T C MB) „ 3 

IT 2 (he) 



(A5) 



where n(e s x) = l/(e e x — 1). Averaging Eqs. (IA3|) and 
over the photon momentum with the distribution 
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function, one finally obtains 



d 2 N (k B T CMB y 2 
■nj{ai)aoii = XQ — arcx 



dtdx 



tt 2 {He) 3 
x / dsPj(s, r y)n(e s x) 



(A6) 



where the redistribution function in Jones's formalism is 
given by 



32 7 4 



32 7 4 



-^- + 2e s + 8 7 2 e 2s 
T 



-4 (A 7 + s) e 5 



87V + 2e 2s 



-4(A 7 -s)e 2s - ^e 3s 



(A7) 



(A8) 



Comparing Eqs. and jMJ with Eqs. (J30J) and (f3T|) . 

respectively, one finds 



P( S)7 )=Pj(s, 7 ). 



(A9) 



which shows the equivalence of the two formalisms for 
extreme-relativistic electrons. 



Appendix B: Derivation of Equations (|57[) and (|58p 

We assume £ ^ 1 for the scattered photons, and 1 <C 
Imin <C for the 7 -parameters. Let us first solve 

Eq. ([3]). It can be rewritten as follows: 



dn(x) 
dr 



2 In 7 max I 'in 

-21n27 mi „ 
2 In 27 mi „ 



—2 In 7 mQS /7„ 

-n(x) . 



n (e 4 7mi „a;j 



dsP IC s, \n [e s 



4 7 ^ 

/7? 



(Bl) 



Then we introduce the following new functions: 

^21nl/_R 

dsP c (s, R) n (e 8 4 7 2 jn x) 

-21n27 mi „ 

= ni(:r) + n 2 {x) , 



(B2) 



21n27 miT , 



21nl/_R 

= n 3 (a;) + n±(x 



dsP IC (s,R)n ( e s — — 



(B3) 



where i? = 7m m/ 7m ax- In Eqs. (|B2I) and (|B3|) . the func- 
tions ni(x), 71.4(0;) are expressed as follows: 



ni{x)= dsP c (s,Q)n(e s Y) , (B4) 

J-21n27 mi „ 
/>oo 

n 2 (x)= dsP c (s,0)n(e s Y) , (B5) 
Jo 

n 3 {x)= f dsP IC (s,0)n(e s X) , (B6) 
J— 00 

<-21n 27 mi „ 

n 4 (x)= dsP IC (s,0)n{e s X) . (B7) 

Jo 

In deriving Eqs. (|B4 [) (|B7|) . we put R = and used new 
variables: 



A 



4 7 2 - ' 

I mm 



(B8) 
(B9) 



Introducing t = e s Y into Eqs. dEJ and |B5]) and /j = e s X 
into Eqs. (|B6j) and (|B7|) . and inserting the explicit forms 
of Eqs. ((40])— ((43]), one obtains as follows: 



ni(x) 



2 7i 2 



3(er-l)- / 

V 'Yj x X a + 5Y' 

1 /ct- 1 , t \ t I 
2 In 



722(2:) 



cr + 3V^ + 3 yyy 

= 6(a-l)(a 2 +4a+ll) ^ (T . , 
(<7 + 1)(ct + 3) 2 (ct + 5) 

dtt-^/ 2 n(t) , 



-ttW), (BIO) 

(T + 1 J 

y(<7-l)/2 

(Bll) 



n 3 (x) 



_ 6(<r- 1)(ct 2 +4ct + 11) 1 1 
" (<r + l)(cr + 3) 2 (o- + 5) ^A^- 1 )/ 2 

x / d^^ 3 )/ 2 ^), (B12) 
Jo 

A 3 r di r 2 



,X d f x dtr 
1 /cr- 1 , t \ t 1 i 2 1 . , . 

Now let us consider the CMB photon distribution func- 
tion 



n a (t) 



1 



e* - 1 



(B14) 



for the initial distribution. Inserting Eq. (|B14|) into 
Eqs. (|BT0j)-(lBT3l), one has for x > 1 



ni(x) = . 
n 2 (x) = . 



(B15) 
(B16) 
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713 (x) 



6(cr- l)(cr 2 + 4cr + 11) 1 1 
(a + 1)(ct + 3) 2 ((7 + 5) 



A' 



1 



Jo e* - 1 ' 

^X 3 f 00 dt 1 r 
ru(a!) = 3(<r - 1)— / — -j — - j-- 



x 3 7^ i e* 

1 /cr - 1 , t \ i 
2 In 



(7 + 3 Vct + 3 



5 

1 i 2 



X X a + lX 2 ) 



(B17) 



1 . (B18) 



The last term in Eq. (IB1|) can be safely neglected for 
x 1. Therefore, one obtains 



dn(x) 
dr 

One finally obtains 
dra(X) _ 1 



713(2;) + 714(3;) • 



(B19) 



dr 



, , f°° dt 1 

*<' -»*•/, 7— 

1 '' < '- 1 -21n.^ ' 



+ 



1 t 2 

n + h ^ (J -\- 3 V cr + 3 " ~ I/I 1 (T + 1I 2 
6(<7 - 1)((T 2 +4(7+11) 1 



(o- + l)(cr + 3) 2 (cr + 5) I( ff - 1 )/ 2 J e *-l 



A - i (cr+3)/2 



(B20) 



One can also solve Eq. ((4]) in a similar manner. It can 
be rewritten as follows: 

dl(x) _ 



dr 



2 111 "irnax /"in 



dsP c s, ) I [e — 

"fmax J \ 47, 
Imin \ T s . 2 



'-2 In 27 mi „ 
r 21n27 mi „ 

dsP/c I s, 

-21n7,„ ax /7 mi „ V Ima. 

-I{x). 



2 

mm 



(B21) 



Then we introduce the following new functions: 

2 In l/R / 

dsP c {s,R)I (e~ s — 

-21n2 7m( „ V 4 7; 

= 7i(a:) + J a (a:), 



2 

mm 



(B22) 



2 In 2 7 „ 



dsP /c (s,i?)/(e- s 4 7 2 



'-21nl/_R 

= h{x)+h{x). (B23) 

In Eqs. (|B22|) and (|B23]l . the functions ^(x), J 4 (a;) 
are expressed as follows: 



2 In 27„ 



dsP c (s, 0) / (p+ s I) 



h{x) = / dsP c (s, 0) / (e^l) , 
Jo 

I 3 (x) = / dsP IC ( S ,0)l(e~ s Y) , 

</ —OO 



(B24) 
(B25) 
(B26) 



h(x) = 



—00 

2 In 27„ 



dsP IC (s,0)l(e- s Y) , (B27) 



where we put R = and used the variables I and Y. 
Introducing t = e~ s I into Eqs. (|B24|) and (|B"2"5j) and 
t = er s Y into Eqs. (|B26|) and (|B27|) . and inserting the 
explicit forms of Eqs. (|4T) |) - (|4"3"|) . one obtains as follows: 



, lW = 3(.-lW'f{-- 



1 fa- 1 



t \ t It 2 



2\ n -) - + ——-L-\l(t)AB28) 

X X a + lX 2 ) y,,y ' 



cr + 3 + 3 

_ , . 6((7-l)((7 2 +4(7+ll) l 

h{x) = 



((7+l)((7 + 3) 2 ((7 + 5) I^-l)/2 
A 



X / dt&^' 2 I{t), 



(B29) 



6((7- l)((7 2 +4(7+ ll) 



3y(a-l)/2 



((7 + l)(a + 3) 2 ((7 + 5) 

/OO 
dtt-^ +7 >/ 2 I(t) , 

* 3 ^dtr 2 i 2 



(B30) 



«*>-«'-i)yj[ £{— 



5 Y 2 



l /ct- l (\ i l 

2 In — l — 



T }/(*). (B3l) 



ff + 3U + 3 Y/Y cr 



Now let us consider the CMB photon distribution func- 
tion 



Jo(t) = /ot 3 n (i) 



(B32) 



for the initial distribution, where lo = (^sTcms) 3 /2-7T 2 
and up(t) is given by Eq. (|Bl4|) . Inserting Eq. (|B32|) into 
Eqs. (IB2H|) - (|B3T|) . one has for x > l 



/ 1 (a;) = 3/o(<T-l)I £ 



/•°° dt l f 

ix- 7e 4 - ll ^" 



l / (7 - 1 , t \ t l f 
2 In ■ 



h{x)=h 



o +3\c7 + 3 ""*Xj X ' cr + II 2 '} 

6((7-l)((7 2 +4(7+ll) l 



(<7 + l)((7 + 3) 2 (a + 5) K— i)/ 2 



A 



X / dt&+^/ 2 - 



l 



e f - l 



h(x)=0, 
h(x) =0. 



, (B33) 



(B34) 



(B35) 
(B36) 



The last term in Eq. (|B2ll) can be safely neglected for 
x ^> 1. Therefore, one obtains 



dl(x) 
~d7~ 



h(x) + I 2 (x). 



(B37) 
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One finally obtains 
dI{X) 



dr 



3(cr 
1 



1)X< 



x 



a - 1 



' dt 1 
Te* - 1 
t \ t 

2 In- 



Comparing Eq. (|B38|) with Eq. (IB20j) . one has the fol- 
lowing relation: 



cr + 5 er + 3\ er + 3 

6(cr - l)(cr 2 + 4a + 11) 1 
(a + l)(a + 3) 2 ( ( 7 + 5) X(--i)/ 2 j 



1 



X t (a+3)/2 

dt— — - -(B38) 



dn(X) _ 1 dl(X) 
dr Iqx 3 dr 

1 1 dI(X) 
" 64 7 ^ m / X3 dr 



(B39) 
(B40) 
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